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We prove that there exists a class of crystalhne insulators, which we call "fragile Mott insulators" 
which are not adiabatically connected to any sort of band insulator provided time-reversal and 
certain point-group symmetries are respected, but which are otherwise unspectacular in that they 
exhibit no topological order nor any form of fractionalized quasiparticles. Different fragile Mott 
insulators are characterized by different nontrivial one-dimensional representations of the crystal 
point group. We illustrate this new type of insulators with two examples: the d-Mott insulator 
discovered in the checkerboard Hubbard model at half-filling and the Affleck-Kennedy-Lieb-Tasaki 
insulator on the square lattice. 



Crystalline insulators are called "band insulators" if 
they can be adiabatically deformed into a system of non- 
interacting electrons in which all the occupied electronic 
bands are separated from empty bands by a finite en- 
ergy gap. Insulating phases with spontaneously broken 
symmetries are clearly distinct phases of matter, which 
cannot evolve into a band insulator without undergoing 
a phase transition. In addition, however, there can exist 
insulating phases of crystalline systems with no broken 
symmetries which non-the-less cannot be adiabatically 
related to a band insulator - such phases are often re- 
ferred to as "Mott insulators" 

Exotic "featureless Mott insulators" or "quantum spin 
liquids" were first proposed Q as candidate ground states 
for highly frustrated quantum antiferromagnetic insula- 
tors. The existence of fractionalized excitations [1, |3| 
and/or topological order (if all the excitations are fully 
gapped) 0, Q are the defining features of such phases. 
For instance, in the 2D RVB state spi'^ 
charge of an electron are separated in the low energy ex- 
citation spectrum. The existence of such phases in vari- 
ous exactly solvable (even though often contrived) mod- 
els has now been established (see e.g. Ref. 0, Q), but 
while there is increasingly interes ting evidence of their 
existence in specific materials IToj . it is fair to say 
that there is still no single material that has been clearly 
proven to possess a spin liquid ground state. 

In this paper, we define a fourth class of crystalline in- 
sulators, which are qualitatively different from the three 
discussed above, which we propose to call "fragile Mott 
insulators." The ground state of such an insulator trans- 
forms non-trivially under the operations of the point 
group. From this it follows that, so long as the allowed 
adiabatic paths respect time reversal symmetry and the 
crystalline point group symmetries, fragile Mott insu- 
lators cannot be adiabatically connected to any band 
insulator; there must be at least one phase transition 
along any path connecting a fragile Mott insulator to a 
band insulator. (For some but not all point groups, it 
is not even necessary to assume time reversal symmetry 



to distinguish fragile Mott insulators [11].) On the other 
hand, fragile Mott insulators are qualitatively different 
from quantum spin liquids in the sense that there are no 
fractionalized excitations and they have a unique ground 
state on the torus. Here we use the term "fragile" in 
describing these phases because they can be sharply dis- 
tinguished from band insulators only when certain crys- 
talline point group symmetries are preserved. We give 
two examples of solvable models of interacting electrons 
which have fragile Mott insulating phases. The protec- 
tion of topological order by crystalline symmetries was 
previously studied (see e.g. Ref. [l^ - [l6| ). For instance, 
its transformation properties under spatial inversion dis- 
tinguishes the Haldane phase of a spin-1 chain from a 
trivial phase isl- 15[. 

To establish that fragile Mott insulators are a distinct 
phase of matter, we rely on the following theorem: 

Theorem: In any time reversal invariant band insu- 
lator, the ground state must transform according to the 
identity representation of the crystal point group. 

Proof: Consider a band insulator described by a non- 
interacting Hamiltonian H where the fermionic opera- 
tor (j>l^ creates a single particle energy eigenstate |0„) = 
(j>l^ |0) . Time reversal symmetry is represented by an anti- 
unitary operator 6, where [Q,'H] = and 9^ — —1, so 
all single- particle eigenstates form Kramers doublets, i.e. 
6(/)„9~^ = (jjfi and 6(/)fi9~^ = -0„, where |0„) and \(f)n) 
are linearly independent states. Since Kramers doublets 
are either both occupied or both unoccupied in a band 
insulator, any band insulator can be described by the 
wave function 
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{n,n)Gocc. 



where "occ." signifies the set of occupied states. Let C g 
G be an element of the point group, G. Since [C, H] — 0, 
any single particle eigenstate can be chosen to be 
a simultaneous eigenstate of C, i.e. C|(/)„) = Xn\4>n)- 
Suppose that the degree of C is m > 1, namely C™ — 1, 
which indicates that its possible eigenvalues are Xji — - 



2 



i2-nj„/m ■ 



1, • • • ,m. Since spatial symmetries com- 
mute with time reversal symmetry, [C, 0] = 0, it follows 

that c |</.ft) = ce \K) = ec \K) = A:e |<^„) = a; |<^^), 

so Xn = A* . With this, we obtain 



c\^o)= n K^fi] i*o> = i*o> 



(2) 



{n,n)Gocc. 



Q.E.D. Note that the above proof is valid for any finite 
lattice size, and remains true as the thermodynamic limit 
is approached. 

The Hubbard Square: The simplest illustration of the 
physics of this new phase comes when we consider the 
ground-state of the four-site Hubbard model with the 
Civ symmetries of a square, for which the most general 
Hamiltonian in the absence of spin-orbit coupling is 



EE{- 



j = l a- 



(3) 



where 



- J 



creates an electron with 



— i^ii:;a,uco an cici^UKjii vviuii SpUl- 

fiolarization a on site j. The ground-state phase diagram 
17, 18] of this "Hubbard molecule" with 4 electrons (in 
average one electron per site) is shown in Fig. [IJb). 

For ^2 > ^1: the ground-state is unique, has spin 0, and 
transforms trivially under the operations of the symme- 
try group, i.e. it has s-wave symmetry. There is a gap to 
the first excited state. At U = 0, the ground-state of the 
non-interacting model is non-degenerate so, in the sense 
of adiabatic continuity, this can be classified as a band 
insulator, and it is so labeled. 

For ti > t2 and U > 0, the ground-state is unique, has 
spin 0, and there is a gap to the first excited state. How- 
ever, the ground-state is odd under rotation by 7r/2 and 
under reflection through various mirror planes - it has 
the symmetry of the (1^2 _y2 representation of the point- 
group. For U = 0, the ground-state is 6-fold degenerate. 



but the structure of the ground-state in the U/ti — ?► 
limit, as we will discuss below, is well defined, and is 
readily seen to be an entangled state that cannot be ex- 
pressed as a single Slater determinant. Because of its 
d-wave symmetry, we call this phase a "d-Mott" state. 

The phase boundary between these two phases occurs 
at ti — t2. At this point, the Hubbard square has a 
higher symmetry - the symmetric group of degree four 
54 or the symmetry of a tetrahedron. Consequently, the 
d-wave and s-wave ground states combine to form a two 
dimensional representation of the symmetric group 6*4; 
the two states are degenerate by symmetry. 

While the Hilbert space in this problem is small enough 
that it is easily solved, it is useful to understand the 
nature of the phases approximately in various limits. In 
the large U/ti limit, the problem is equivalent to the 
spin- 1/2 Heisenberg square, with exchange coupling Ji — 
A£\/U and J2 = '^t^/U between nearest and next-nearest 
neighbor spins, respectively. There are two spin-0 states. 
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where 6-^- = [c^..^cj^ + cj^c^^]/-\/2 creates a singlet pair 
between sites i and j. Each of these states is an eigenstate 
of the Hamiltonian, since \^ s) is the unique spin-0 state 
which transforms according to the identity representation 
of the point group, while \^d) is the unique state with 
dj.2_y2 symmetry. While \^s) is a simple valence bond 
state, \'^d) is a "short-ranged resonating valence bond 
state," along the lines that were considered in Ref. il,[5j. 

Still more interesting is the small U limit. For U = Q 
and t2 > fi, the ground-state is readily seen to be the 
simple Slater determinant state. 
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the 



ground state in the U /ti 



j ^. However, for ti > t2, 
limit is the entangled state 
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FIG. 1: (a) Schematic representation of the checkerboard lat- 
tice. Intraplaquette hoppings are ti between nearest neigh- 
bors and t2 along diagonal direction and interplaquette t' . (b) 
Phase diagram of the Hubbard square (solid lines) and of the 
checkerboard Hubbard model (dashed lines) in the small t' 
limit. 
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The fact that even in the U/ti — > limit, the ground- 
state is entangled is a remarkable feature of this partic- 
ular Hubbard molecule which is directly related to 
the existence of the d-Mott phase. 

The d-Mott insulator. The existence of a d-Mott phase 
in an extended system follows simply from the solution 
of the Hubbard square [l^ . Consider the "checkerboard 
Hubbard model," shown in Fig.[TJa), which we can think 
of as a model of a molecular crystal built of Hubbard 



squares 
j 



Now, c 



t 

aj,cr 



creates an electron on the 



site of the a*'' molecule with spin polarization cr, and 



{o.j,a'j' 
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where H^'^ is the molecular Hamiltonian defined as in 
Eq. with Cj ^ — > ^- ^ and {aj, a'j') signifies pairs 
of nearest-neighbor sites on adjacent molecules. In the 
t' — ?► limit, the system consists of disconnected square 
molecules so the ground-state is given by the tensor prod- 
uct of the ground-states of the isolated squares. It fol- 
lows from the existence of a gap in the spectrum that the 
thermodynamic properties of the system evolve smoothly 
for non-zero t' as long as t' is small enough and can be 
computed perturbatively. In particular, all symmetry 
properties of the ground state are independent of t' in 
this regime, so with the possible exception of a region 
of width ~ t' /ti about the line ti = t2 (which we dis- 
cuss separately, below), the phase diagram in Fig. ^h) 
applies to the thermodynamic limit of the checkerboard 
Hubbard model for small enough but non-zero t' /ti. 

It remains to establish, for the extended system, 
whether or not each phase can be adiabatically connected 
to a non-interacting insulator. The band-structure of this 
model with t/ = is easily computed. For t2 — ti > \t'\, 
there is a band-gap separating the occupied and empty 
states. As this state is approached smoothly as C/ — 0, 
it is still clearly the case that the label "band-insulator" 
applies to the large ^2/^1 phase. 

To prove that the d-Mott phase at ^1 — ^2 > \t'\ is not a 
band-insulator, we consider a system of L by L molecules 
with periodic boundary conditions and L odd. Then, it 
is obvious that the ground state wave function acquires 
a minus sign under 90° rotation (C4) around a plaquette 
center or refiections (labeled by ad) along a diagonal line 
across the plaquette centers. 



l^'d-Mott) = — l^'d-Mott) ,C4 l^'d-Mott) 



I* 



d-Mott 



,(8) 



which implies the ground state wave function transforms 



like a 



orbital under the point group C^y . Accord- 



ing to our theorem, this implies that the d-Mott phase 
cannot be adiabatically related to a band insulator. If 
L is even, the ground-state transforms according to the 
trivial representation of the point group. Thus, for these 
boundary conditions, we cannot use symmetry to prove 
that the system is not adiabatically related to a band in- 
sulator. However, the identity of a phase should not de- 
pend on the way the thermodynamic limit is approached. 
Thus, we believe that proving that adiabatic evolution to 
a band insulating state is impossible for L odd is suffi- 
cient for establishing that this is a distinct phase (23| . 

To gain further insight into this problem, we focus on 
the behavior of the system near the boundary between 
the band and d-Mott insulating phases. Here, for t' = 
and ti = t2, there are two low- lying states per molecule, 
so to study the effects of non-zero coupling t' between 
molecules and non-zero \ti — ^2!, we must derive an ef- 
fective Hamiltonian using near-degenerate perturbation 
theory. We thus define pseudo-spin operators associated 
with each molecule, where rf^ = 1 if molecule a is in 



its s-wave ground-state, and = — 1 if it is in its d- 
wave ground state. We define the raising operator r+ 
that promotes the molecule from its d-wave to its s-wave 
ground state, etc. In terms of these operators, we de- 
rive a pseudo-spin i quantum Ising model defined on the 
molecular lattice to lowest order in 1^1 — ^2! and to second 
order in t': 



J. 



eff 



(Ta ■ e+){fa+x ■ e+) 



+(ra-e_)(' 



a+y 



(9) 



where e± = -^i;±iz are unit vectors, Jeff = j^gi{U/ti), 

and h = {t2 - ti)g2{U/ti) + ^gsiU/ti) where 31,2,3 
are regular non-negative functions of U/ti which can 
be computed numerically. Under a C4 rotation of the 
system, the lattice points transform as a — >■ a', and 
Ta — > e^'^'^'^/^Ta'e""'^^/^; the Hamiltonian is invariant un- 
der this transformation. Any ordered state with a non- 
zero expectation value of breaks the global C4 sym- 
metry of the lattice. For h = 0, the model has two ad- 
ditional global Ising symmetries: Ta gi-^T^ /^f^^—nrT^ /2 
and Ta — > e^'^'^^'/^fa'C""^'^"'/^. For h > J^s, the ground 
state is polarized with t^ = 1 (band insulator) while for 
h <^ — Jeff, the ground state has t^ — —1 (d-Mott). We 
do not, as yet, have definitive results on the nature of 
the transition(s) between these two extremal phases for 
small \h\. Based on the obvious similarity between Hcs 
and the transverse field Ising antiferromagnet, it seems 
likely that there is at least one intermediate phase with 
spontaneously broken translational and C4 symmetry oc- 
curs in a region of width ~ i'^/ii about the ti = t2 line 
as shown by the dashed lines in Fig. [TJb). 

The AKLT insulator. As a second example, we will 
consider a model with a low energy effective theory equiv- 
alent to the famous AKLT model of quantum antifcrro- 
magnetism [24]. For simplicity, wc define the model on 
the same checkerboard lattice as shown in Fig. la with a 
mean electron density of one electron per site (4-electrons 
per unit cell), but with the Hamiltonian 

H = A-^(6 - + J E ^4(a, a') + A ^(5,^)^(10) 

a {a.a') a 

where the total spin on square a is 



'3 = 1 ^a..j 



c') and P4(a,a') is the pro- 



jection operator onto spin-(Sa -|- Sa')^ = S{S -\- 1) with 
S — A. This Hamiltonian has all the same spatial sym- 
metries as the checkerboard lattice. 

For if, J > 0, and A = 0, this Hamiltonian is positive 
semi-definite. The K term is minimized by any state 
with 4 electrons in a spin 2 state on each molecule. In the 
ground-state subspace defined by this K term, the J term 
in Eq. ([T0| , is precisely the spin 2 AKLT Hamiltonian on 
the square lattice 2J|. The zero energy ground state is 
thus seen to be 
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KLT/ — 



E 



1 Jo,2 Ja,3 



|0) 



,(11) 



where (ja,i, Ja,2, Ja,3, ia,4) labels one of the 24 possible 
permutations of (1, 2, 3, 4) on each plaquette a and (—1)''" 



is the signature of permutation. Here, 



I a/2 creates 



a singlet pair between two sites (a,ja,i) and (a+a;, 
and fet „ „ , - „ is defined similarly. A/" is a normal- 
ization constant. Again, for the system on a L by i 
lattice with L odd, we obtain 

0-i,(d) I^AKLt) = -|*AKLT>, (12) 

which indicates that the AKLT insulator on the checker- 
board lattice transforms as the A2-representation under 
the C4t, group and thus cannot be adiabatically con- 
nected to any band insulator. Moreover, since the d-Mott 
insulator transforms in the i?i-representation (dx2-y2), 
as discussed previously, the AKLT insulator and d-Mott 
insulator are distinct states of matter. Different fragile 
Mott insulators are distinguished by the different one di- 
mensional representations of the point group in question. 

For K, A > 0, and J = 0, the system becomes dis- 
connected squares, each of which has a unique ground 
state with S — 2 and — 0, and the system's ground 
state is then given by the direct product state |2?) = 
Yla ^\Sa — 2, — 0). Due to the existence of a gap in 
the spectrum, the qualitative nature of I'D) survives for 
non-zero but small enough J. Since I'D) on a L by L lat- 
tice with L odd transforms like the d2;2_j^2-representation 
under Ciy, it is adiabatically connected to the d-Mott 
state but not to the AKLT insulator. 

Concluding remarks: Beyond the point of principle - 
that fragile Mott insulators are thermodynamically dis- 
tinct from band insulators - there is little to distinguish 
them in practice. In terms of their elementary excita- 
tions, they are indistinguishable from a simple semicon- 
ductor. A proposal was made in Ref. [3l for a "phase- 
sensitive" measurement of a d-Mott phase, involving use 
of a fragile Mott insulator as the insulating component 
of superconductor-insulator-superconductor junctions of 
various geometries. While the experiments involved may 
not be simple, the analysis highlights another perspective 
on this new phase. 

Let us again consider the checkerboard lattice, initially 
in the limit t' — >■ 0. We define a new vacuum state to be 
the tensor product of the ground-states of each square 
with two electrons per square. We can therefore view 
17j the ground-state of the extended system as being a 
Bose insulator with a single d-wave pair per square. The 
same Bose insulating state made of s-wave Cooper pairs 
would likely be adiabatically connected to a band insu- 
lator, but in the d-wave case, the insulating state retains 



information about the internal structure of the pairs. In 
the above mentioned experiments, this gives the same 
sign of any phase-sensitive measurement that would be 
obtained using a piece of d-wave superconductor as a link 
in a SQUID loop although of course in the case of the in- 
sulator, the Josephson coupling is strongly attenuated. 
Needless to say, as with all difficult to realize theoretical 
proposals these days, a d-Mott phase could be realized 
for cold atoms in optical lattices under appropriate cir- 
cumstances [iil. l26l|. 
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